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Taste breaking effects in scalar meson correlators 
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As a consistency check of the staggered-fermion fourth-root approximation, we analyze the arj 
and fo correlators, including the effects of two-meson taste-multiplet intermediate states. Rooted 
staggered chiral perturbation theory describes the contributions from the pseudoscalar taste mul- 
tiplets in terms of only a few low energy constants, which have all been previously determined 
by the MILC collaboration. In previous work one of us (Prelovsek) showed that the two-meson 
"bubble" contributions could explain the observed anomalies in the lattice data for the isovector 
flo channel. In the present work we extend this analysis to the fy channel. On a MILC collab- 
oration lattice ensemble at 0.12 fm with 2+1 flavors of Asqtad-improved staggered fermions, 
we have made new measurements of correlators in both channels for a variety of momenta. A 
fit to these correlators gives low energy constants that are reasonably consistent with previous 
determinations by the MILC collaboration. 
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1. Introduction 

When up and down quark masses are sufficiently light, the isovector (oq) and isosinglet (/o) 
scalar meson correlators are dominated at large distances by two-body states composed of n, K, 
and v\ . In the lattice staggered fermion formulation, residual fermion doubling results in a sixteen- 
member taste multiplet for each meson. When this symmetry is broken at nonzero lattice spacing, 
the multiplet is split, resulting in a proliferation of two-body intermediate states, thereby complicat- 
ing the analysis of the correlators. To make matters worse, many of these states are lattice artifacts 
with unphysical masses, and some have ghost (negative) weights. Thus we see nonlocalities in the 
form of unphysical long-range contributions to the correlators, and we see violations of unitarity 
in the form of unphysical and negative norm intermediate states. Such artifacts are expected to 
disappear in the continuum limit. 

Fortunately, these complexities are described by staggered chiral perturbation theory (S^PT). 
The rooted version of this theory (rS^PT) ^ ||] provides a strict framework in which to analyze 
these two-body contributions in terms of just a few low energy constants - constants that are also 
determined in an independent staggered chiral analysis of the pseudoscalar meson masses and 
decay constants [J3j] . 

Thus the study of scalar correlators provides an explicit example of the anticipated diseases 
of the fourth-root approximation and the manner in which they are cured in the continuum limit. 
Further, it provides another test of the ability of rS^PT to model rooted staggered fermion QCD. 
Previous tests include the aforementioned fits to the pseudoscalar meson masses and decay con- 
stants and a comparison of the topological susceptibility measured in rooted staggered fermion 
QCD and in rS^PT[g||. 

2. Pseudoscalar meson multiplet 

The cto channel was studied in recent years in staggered fermion QCD by the MILC collabora- 
tion [j|] and UKQCD collaboration [EJ. Both groups found that the correlator appeared to contain 
states with energies well below possible combinations of known mesons. At Lattice 2005 one of 
us showed that rS^PT provides a simple explanation [0 §J. In Fig. [j] we show taste splittings of 
the low-lying pseudoscalar mesons for the MILC a = 0.125 fm ensemble, as predicted in rS^PT 
In the case of the pion the predictions were confirmed in lattice measurements [§, ^]. Of particular 
concern are the tj and tj' multiplets. They are mixed by the anomaly. But the anomaly is a taste 
singlet. So only the taste-singlet members are proper candidates for the physical states. The other 
taste members are unmixed at tree level. We have arbitrarily labeled the lighter member tj and 
heavier member tj'. The taste-axial- vector and vector hairpin couplings mix the 77 's weakly. We 
ignore that effect in the figure, but take it fully into account in our analysis. 

Now, let us examine the various two-meson contributions to the ao correlator. In proper QCD 
the lowest intermediate state has mass m n + m-q. In staggered QCD the physical ao is a taste 
singlet, so we need to examine the various two-meson contributions to the taste-singlet correlator. 
Taste symmetry requires that both pseudoscalar mesons in an intermediate state have the same 
taste assignment. In particular, the nrj intermediate states include a taste-pseudoscalar pion and a 
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Figure 1: Taste multiplets for K, K, rj, and r\' for the study ensemble. The letters indicate the taste symmetry 
assignments, grouped as predicted by Lee and Sharpe [|To||. 

taste-pseudoscalar eta with a mass equal to twice the Goldstone pion mass. Other taste pairs are 
similarly light. Thus we expect an anomalous low energy contribution to the correlator. 

3. Threshold weights 

Explicit expressions were presented in [IvL [S|] for the two-meson contribution to the ao correla- 



tor, including all tastes. We have extended this analysis to the f correlator [11]. (See also Sec. VI 
of Ref. [^] for a description of the result in the one-flavor case.) The expressions (not written here) 
consist of several two-body threshold contributions with weights depending on the number of fla- 



vors and the replica factor [12] needed to correct for the unwanted taste multiplicity. If we ignore 
mixings induced by the taste axial-vector and vector hairpins, the weights for the physical 2 + 1 
flavor case are given in Table [I]. In the continuum limit in the oq channel, all thresholds but the 
taste-singlet 7ir\ become degenerate and the axial-vector and vector hairpin mixings vanish. As is 
evident from the table, the weights add up to zero. The ghost threshold has eaten the unphysical 
thresholds, and only the physical taste-singlet %r\ survives. In the fo channel, all the %% thresholds 
become degenerate leaving a proper weight of 3 for the three charge assignments for the physical 
pions. 

For the single-flavor case there is no Goldstone boson, since the one pseudoscalar is lifted by 
the anomaly. Yet at nonzero lattice spacing there is a multiplet of near Goldstone bosons composed 
of taste components that are unmixed by the anomaly. In the continuum limit the two-body thresh- 
olds for these states become degenerate and the weights add to zero in the fo channel as can be 
seen from Table ^[ Thus the unphysical states decouple from this channel. 

4. Fits and Results 

Measurements were done on a single ensemble of 510 gauge configurations of dimension 
24 3 x 64, generated by the MILC collaboration with 2+1 flavors of staggered quarks of bare mass 
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Table 1: Two-body threshold weights for 2+1 flavors in the «o an d fo channels. The «o weights apply to 
the %r\ contributions for tastes I,P,V,A, T, and the fo weights apply to the %% contributions. Here the taste- 
singlet r\ is the physical state. For the ao the weight 10 denotes the contribution from the taste-singlet pion 
plus the bare taste-singlet rj (unshifted by the anomaly). For the fo the weight 10 denotes the contribution 
from two bare taste-singlet etas. 
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Table 2: Two-body threshold weights for the single flavor case fo channel. The notation is the same as Table 
0- 

taste fo 

"To -15/8 

V 4/8 

T 6/8 

A 4/8 

P 1/8 



amud = 0.005 and am s = 0.05 at a lattice spacing of approximately a = 0. 125 fm. For this ensemble 
the pi/rho mass ratio is 0.303. 

We measured the correlator of scalar densities ij/y for the isovector and isoscalar channels and 
for five low-lying total momenta p. The correlators were fit to the expressions 

C a o{pf) = fmeson,ao{Pf)+fbubble,ao{pf) 
CfoiPf) = f meson, fo(pf) + /bubble, fo(pf) 

and where the single-meson contribution is modeled by a single pole of each parity, 

fmesonMpf) = b a0 (p) exp [-E a0 (j>)t] + b nA (/?)(-)' exp[-£ M [p)t] + (N t -t) 
/meson,/o(/V) = c 8(p, 0) + b f0 {p) exp[-E f0 (p)t] + b 1]A {p){-) t exp[-E ntA (p)t] + (N t -t), 

and where the "bubble" or two-meson contributions are dictated by rS#PT[g, 0, |llj|. The meson 
contributions require 13 parameters (four for the masses, two each for the momentum dependence 
of the four coefficients, and one for the zero-momentum vacuum disconnected term in the fo). For 
each meson energy we assume a standard dispersion relation E(p) 2 = m 2 + p 2 . The meson masses 
in the two-meson bubble contributions are fixed from other spectroscopic studies plus the known 
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taste splittings ||5|,[3j]. There remain three low energy constants from rS^PT, namely, 

M =m n /{2m U 4), 8 A =a A 8' A , 8 v = a A 8[ / 

All three parameters were previously obtained from a study of meson masses and decays |^. We 
found it necessary to fix (through a Bayesian prior) the vector hairpin constant 8y to the value 
obtained in that study. The remaining two constants were then adjusted with only loose priors. 
The resulting best fits are shown in Figs. ||, ||, and ^. 

Table 3: Comparison of our fit parameters for the rS^PT low energy constants with results from [||] 

quantity our result meson mass, decay 

n m 2 n /(2m u4 ) 8.2(1.1) 6\7~ 

5y (prior) -0.016(23) 

8 A -0.056(7) -0.040(6) 

The agreement is worse if we used the value of n/i in the three-flavor chiral limit (« 4.5) 
from those fits, rather than the higher-order m|/(2m M ^), suggesting, perhaps, that a higher order 
calculation of the bubble contribution might improve the agreement. 
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Figure 2: Fit to the ao correlator for five momenta. 

5. Conclusion 

We have found that rooted staggered chiral perturbation theory together with a couple of ex- 
plicit low-lying mesons provides a reasonably consistent accounting of the ao and fo correlators 
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Figure 3: Fit to the /o correlator for four momenta. 
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Figure 4: Fit to the /o correlator for zero momentum. 
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for several low momenta. The spurious long-range spectral components are lattice artifacts, which 
the theory predicts will disappear in the continuum limit. 
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knowledged. This work was supported by the US NSF and DOE. We thank the MILC collaboration 
for the use of its lattice ensembles. 
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